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We study accelerating relativistic reference frames in Minkowski space-time under the harmonic 
gauge. It is well-known that the harmonic gauge imposes constraints on the components of the metric 
tensor and also on the functional form of admissible coordinate transformations. These two sets of 
constraints are equivalent and represent the dual nature of the harmonic gauge. We explore this 
duality and show that the harmonic gauge allows presenting an accelerated metric in an elegant form 
that depends only on two harmonic potentials. It also allows reconstruction of the spatial structure 
of the post-Galilean coordinate transformation functions relating inertial and accelerating frames. 
The remaining temporal dependence of these functions together with corresponding equations of 
^ | motion are determined from dynamical conditions, obtained by constructing the relativistic proper 

^vq . reference frame of an accelerated test particle. In this frame, the effect of external forces acting 

on the observer is balanced by the fictitious frame-reaction force that is needed to keep the test 
' particle at rest with respect to the frame, conserving its relativistic linear momentum. We find that 

(^S , this approach is sufficient to determine all the terms of the coordinate transformation. The same 

method is then used to develop the inverse transformations. The resulting post-Galilean coordinate 
00 ' transformations extend the Poincare group on the case of accelerating observers. We present and 

discuss the resulting coordinate transformations, relativistic equations of motion, and the structure 
of the metric tensors corresponding to the relativistic reference frames involved. 
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I. INTRODUCTION 



Most modern precision physics experiments are conducted in non-inertial reference frames corresponding to either 
the surface of the rotating Earth or an accelerating spacecraft [l|, Q . Quite often, in addition to special precautions 
aimed at reducing the ambient non-gravitational acceleration noise, one also needs to develop a relativistic treatment 
of all the observable quantities. Unfortunately, there is no general agreement on how to approach these problems in 
a relativistic manner. For instance, the treatment advocated in [3|— 15| is incompatible with most of the methods of 
relativistic reference frames [Ml and also is inadequate from the practical standpoint involving the transformation 
between the experimental (accelerated) and laboratory (inertial) frames. In general, it is not clear how to write a 
suitable explicit form of the metric tensor of the accelerated reference frame and construct coordinate transformations 
between the inertial and accelerated frames, especially when high accuracy is required. 

Within the realm of special theory of relativity, to describe the dynamics, one uses the relativistic mechanics of 
■ Poincare which uses the Lorentz transformations to describe physical process in various inertial reference frames: 

c_ 7 ( (+ i^i), ✓_, + ^ +(T _ 1) i8<ip), 7=(i-!)~ 4 . a) 

Note that the transformations inverse to those of ([TJ arc obtained by simply replacing (t' , r') — s- (t, r) and the velocity 
v — > — vq. The Lorentz transformations are well suited to study the case of uniform motion with constant velocity. 
However, while analyzing experimental data one often has to deal with acceleration and dynamical noise whose 
presence limits the practical applicability of these results. To address the presence of acceleration, for each instant 
of time one can use (UJ to conceptually define a set of instantaneous quasi-inertial reference frames. Although useful 
conceptually, this approach is insufficient in practice as it leads to a decreased precision on larger time scales. 

When the Riemannian geometry of general theory of relativity is concerned, it is well known that coordinate charts 
are merely labels. Usually, space-time coordinates have no direct physical meaning and it is essential to construct 
the observables as coordinate-independent quantities. Thus, in order to interpret the results of observations or 
experiments, one picks a specific coordinate system, chosen for the sake of convenience and calculational expediency, 
formulates a coordinate picture of the measurement procedure, and then one derives the observable out of it. It is 
also known that an ill-defined reference frame may lead to appearance of non-physical terms that may significantly 
complicate the interpretation of the data collected @ . Therefore, in practical problems involving relativistic reference 
frames, choosing the right coordinate system with clearly understood properties is of paramount importance, even as 
one recognizes that in principle, all (non-degenerate) coordinate systems are created equal [l5j . 
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Modern theories of rclativistic reference frames [6l-ll3| , dealing predominantly with the general theory of relativity, 
usually take the following approach: As a rule, before solving gravitational field equations, four restrictions (coordinate 
or gauge conditions) are imposed on the components of the Ricmannian 1 metric g mn . These conditions extract a 
particular subset from an infinite set of space-time coordinates. Within this subset, the coordinates are linked by 
smooth differentiable transformations that do not change the coordinate conditions that were imposed. A set of 
differential coordinate conditions used in leading theories of relativistic reference systems, such as that recommended 
by the International Astronomical Union (see, for instance [HI), are the harmonic gauge conditions. The harmonic 
gauge has a very prominent role in gravitational physics, starting with the work of Fock [l6|. In addition, a set of 
specific conditions designed to fix a particular reference frame is added to eliminate most of the remaining degrees of 
freedom, yielding an explicit form for the coordinate system associated with either inertial or accelerating frames. 

In the present paper we introduce a new method for deriving the metric associated with the proper reference 
frame of an accelerated observer. To constrain the available degrees of freedom, we will also use the harmonic gauge. 
An argument in favor of choosing the harmonic gauge is that tremendous work in general relativity has been done 
with this gauge, which was found to be a simplifying and useful gauge for many kinds applications [Hj]. Using the 
harmonic gauge, we develop the structure of both the direct and inverse coordinate transformations between inertial 
and accelerated reference frames. Such mutual transformations are not fully treated in the current literature, since 
usually either the direct @, HH or the inverse [l2[ transformation is developed, but not both at the same time. The 
method presented in this paper could help to generalize contemporary theories of relativistic reference frames and 
to allow one to deal naturally with both transformations and relevant equations of motion in a unified formalism. 
Finally, the method presented here does not rely on a particular theory of gravitation; instead, it uses the covariant 
coordinate transformations to explore the dynamics in Minkowski space-time from a general perspective. Thus, any 
description from the standpoint of a metric theory of gravity must yield our results in the special relativistic limit. 

We begin our discussion in Section [TT] by generalizing the Lorentz transformations (TTJ) to the case of accelerated 
motion. This generalization can be done in the form of a slow motion approximation that admits expansion of all 
the quantities involved in the form of power series. This task is accomplished by introducing acceleration-dependent 
terms in the coordinate transformations. The resulting transformation is given in a general form that relies on a set 
of functions that are precisely determined in the subsequent sections. 

The local coordinate system of an accelerated observer is not unique. We use the harmonic gauge to constrain 
the set of coordinates chosen in the accelerated reference frame. We carry out this task in Section [IIII and determine 
the metric tensor describing the accelerated reference frame and the structure of the coordinate transformations that 
satisfy the harmonic gauge conditions. We observe that the metric in the accelerated frame has an elegant form that 
depends only on two harmonic potentials, which yield a powerful tool that allows reconstruction of the spatial part 
of the structure of the post-Galilean coordinate transformation functions between inertial and accelerating frames. 

To fix the remaining degrees of freedom and to specify the proper reference frame of an accelerated observer, 
we introduce a set of dynamical conditions in Section IIVI Specifically, we require that the relativistic linear three- 
momentum of the accelerated observer in its proper reference frame to be conserved. This conservation leads to fixing 
the time-like coordinate in the accelerating frame, which allows to fix uniquely all the remaining terms. 

A similar approach can also be carried out in reverse, establishing coordinate transformation rules from an acceler- 
ating to an inertial reference frame, which is done in Scc.fV] Key to the approach presented in this section is the use of 
the contravariant metric for the accelerating frame, which leads in a straightforward manner to the inverse Jacobian 
matrix. This allows us to present the inverse transformations, in which the roles of the inertial and accelerating 
coordinates are reversed. The calculations are formally very similar to those presented in Sees. IIIII and HVl thus, to 
avoid repetitiveness, we show only the main results. 

We conclude by discussing these results and presenting our recommendations for future research in Section fVTl We 
also show the correspondence between our results and those obtained previously by other authors. 

II. POST-GALILEAN COORDINATE TRANSFORMATIONS AND BOOKKEEPING 

To describe the dynamics of the TV-body problem one usually introduces TV + 1 reference frames with their own 
coordinate charts. We need one global coordinate chart defined for the inertial reference frame that covers the entire 



1 The notational conventions employed here are those used by Landau and Lifshitz |14H : Letters from the second half of the Latin alphabet, 
m, n, ... = 0...3 denote space-time indices. Greek letters a, /?,... = 1...3 denote spatial indices. The metric 7 mn is that of Minkowski 
space-time with 7, n n = diag(+l, — 1, — 1, — 1) in the Cartesian representation. The coordinates are formed such that (ct, r) = (x°,x a ), 
where c is the velocity of light. We employ the Einstein summation convention with indices being lowered or raised using ■jmn- 
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system under consideration. In the immediate vicinity of each of the N bodies in the system we can also introduce a 
set of local coordinates defined in the frame associated with this body. In the remainder of this paper, we use {a;™} to 
represent the coordinates of the global inertial frame and {y m } to be the local coordinates of the accelerated frame. 

The most general post-Galilean form of a finite transformation between the coordinates of inertial {x m } and 
accelerating {y m } dynamically non-rotating reference frames, may be given in the following form [l7j : 



x 



?y° + c~ 2 JC{y\ f) + c- 4 £(y°, + 0(c~ 6 ), (2) 
= y» + zZ(y )+c- 2 Q»(y ,yt) + O(c- 4 ), (3) 

where Zq is the Galilean vector connecting the spatial origins of two dynamically non-rotating frames, and we introduce 
the post-Galilean vector Xq(j/°) connecting the spatial origins of the two frames 

a#(y°) = % + c- 2 Q»(y°, 0) + O^). (4) 

The functions K, C and Q M are yet unknown. It is anticipated that these functions depend only on the relative motion 
between the reference frames involved. Our objective is to construct an explicit form of the proper reference frame of 
the accelerating observer by determining the explicit functional form of the coordinate transformation functions K., C 
and Q M . We shall refer to this approach as the K.CQ- formalism for rclativistic reference frames. 

In this formulation, we use the dimensioned parameter c~ 4 as a bookkeeping device for order terms. For instance, 
when we write c~ 2 K. in Eq. @ above, this implies that the function K is of order v 2 y° with v -C c being the 
characteristic velocity of an observer; hence, c~ 2 JC is of order (v/c) 2 y°, which remains small relative to y°. Similarly, 
c~ 4 £ ~ (v/c) 4 y° and c _2 Q M ~ (u/c) 2 ^. Specifically, in the no-acceleration limit, Eqs. (JS])^© should reduce to the 
Lorentz transformations given by Eqs. (JT|), which can be written in approximate form as: 



x = y+^ + ^jvoc-V + ^^ + Ofc- 4 ), (6) 

corresponding to the following set of ICC Q coordinate transformation functions: /Co = 1 ;U 2 ?/ -l-c(vo-y)-|-0(c _4 ), Co = 
|u 4 ?/ + |u 2 c(v • y) + 0(cT 2 ), 0% = ±v (v • y) + §v§v cr V + 0(c" 2 ), with = Voc" V- These approximations 
remain valid so long as y^/y° <C vq/c, a condition that is naturally satisfied along the world-line y^ = 0. 

The coordinate transformation rules for the general coordinate transformations above are easy to obtain and express 
in the form of the Jacobian matrix dx m /dy n . From Eqs. (JSJ-Q, we get: 

dx° _ odfC A dC fi , dx° , dlC A dC „. 

w = 1 + c w + c W + ° {c } ' w = c w + c w + ° {c 1 (7) 

where Vq = Zq = cdzy/dy is the time-dependent relative velocity between the two frames. 

The transformations given by Eqs. ([I])-© must be smooth and the corresponding Jacobian matrix dx m /dy n must 
be non-singular, i.e., its determinant must be nonvanishing 2 : 

**(w) - ,+ '- 1 (3f + «« v '° + w) + ° <<r4) * <»> 

which guarantees the existence of inverse transformations, discussed in Sec. |V1 



2 We can find an upper limit for the distances at which the condition remains valid. With the help of Eqs. (1 3 2 1 ) — (T33T l and the solution 
for the function K given by (1971) below, we can present Eq. ifSJ as 

det (^r) = 1 - ^ + M) + ^~ 4 ) = 1 + ~, (W) ~ <*>) + °( c " 4 )- 

Therefore, the transformations l(2|-|(3]l arc non-singular up to the distances from the world-line that satisfy the following condition 
y < (c 2 — 2ipo)/(2ao). Note that, for any realistic scenario, this condition is not violated on solar system scales 
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To establish the form of the metric tensor in the moving frame, we apply the usual tensor transformation rule: 

dx^ dx^ 

Vmn(y) = QymQyn ^Mv)), (10) 

which, together with Eqs. (J7J) (JSJ) ? allows us to compute the metric components in the arbitrarily moving frame: 

-, 2 (dlC , f i 2 ( dC dQ e a ,fdJC\^^ fi , 

voo(y) = 1 + ^{^ + H^} + ^{^ + -o e ^ + i(^ T ) } + 0(c-% (11) 

lrl 3K i lfl3£ 9Q £ dQ x 1 die die ^ _ B , 

I (1 dlC 1 dlC dQ x 8Q X \ _ 4 . 

Eqs. (|llj) -([13 p show the structure of the metric tensor of the Minkowski space-time in an arbitrarily moving reference 
frame. The actual dependence of rj m n on the transformation functions (/C, L, Q a ) will be important to define the proper 
relativistic reference frame of an accelerated observer under the harmonic gauge. This will be investigated next. 

III. IMPOSING THE HARMONIC GAUGE CONDITIONS 

The dynamical condition, i.e., the requirement that the spatial origin of the transformed system of coordinates 
x m _ J(y m ) i s to move along a specific world-line, does not uniquely define y m . The existence of this coordinate 
freedom allows us to impose the harmonic (de Donder) gauge condition in the local frame 3 : 

drn(V=VV mn ) =0. (14) 

The vanishing of the covariant derivative of the metric tensor allows us to present Eq. (|14j) in the following equivalent 
form: 77^1^(77) = 0. Remembering the transformation rules of the Christoffel symbols under general coordinate 
transformations [l4[ , one can verify that the latter equation is equivalent to imposing the harmonic conditions on the 
transformation functions in Eqs. <(2j)— (J3j) : 

U v x m = 0, (15) 

where D r) = {^—i])^ 1 d m (^—rji] mn d n ) denotes the covariant d'Alcmbcrtian, in this case acting on x m , which are 
treated as individual scalar functions. 

Therefore, on the one hand the harmonic gauge imposes restrictions on the partial derivatives of the metric tensor, 
as seen in Eq. (fT4l . On the other hand, it restricts the choice of admissible coordinate transformations only to those 
that satisfy the harmonic equation (|15[) . These two consequences of imposing the harmonic gauge will be used to 
establish the structure that the metric tensor (|11[) - (|13[) must satisfy under the harmonic coordinate transformations 
©-(131) and also to constrain the form of the transformation functions (/C,£, Q a ). 

A. The form of the metric tensor in a moving frame 

We require the proper reference frame of a moving observer to exhibit no rotation or shear of its coordinate axes. 
The components of the metric tensor r\ mn that satisfy these requirements under the harmonic gauge conditions Q14[) 
must, therefore, satisfy the following set of partial differential equations 

i£j = G(c- 4 ), (16) 
\cd {^ - 7 £ \S } + d v r$l = 0(c- 2 ), (17) 



3 In Sees. Hill and HVl when convenient and unambiguous, we shall use the abbreviated form of the partial derivative operator to represent 
partial derivatives with respect to {y m } : d m = d/dy m . 
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By formally integrating Eq. (|T8)) and choosing the solution corresponding to a subset of the harmonic gauge with 
spatially isotropic coordinates, we are led to the following form of the gauge conditions ([T5)) - ([T5| : 

r)W = 0(c- 4 ), (19) 
2cd ^+d^ = 0(c~ 2 ), (20) 
*?S+Woo = 0{c- 2 ). (21) 

Eqs. ([T9"|) - (|2"Tj) represent a set of harmonic conditions on the metric tensor r) mn given by Eqs. (fTT |) — (|T5|) in the coor- 
dinates of a moving reference frame. This set of gauge conditions forms the foundation of our method of constructing 
a proper reference frame of an arbitrarily moving observer. As a result, the metric representing space-time in the 
moving frame may be presented in the following elegant isotropic form that depends only on two harmonic potentials: 

MV) = l-^u(y) + ^u 2 (y) + 0(c~ 6 ), (22) 
Vo a (y) = - la Au\y) + 0(c' 5 ), (23) 
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Vap(y) = J a i3+l a /3-^u(y) + 0(c ), (24) 



where we introduced the scalar u(y) and the vector u a {y) potentials, defined as: 

i„„.„? } k,~7!n, 1- — I — 



/ \ i e 1 rot cy ifoi^v (ok i ,y\ _ 4 , 

= 2 u o^o-^(^o +"'0^+ a(j£oJ - (^0 + 2^0 J j + 0(c ), (25) 



5Q e 




le dy° 


2 \dy°J 






dy° 



It follows from Eq. ((20)) that these potentials satisfy the continuity equation: 

cdou + d.u' = 0{c- 2 ). (27) 

With the help of Eqs. (|2"9")) - (P3")) . derived in Sec. MI Bl below, we can verify that the potentials u and u a satisfy the 
following harmonic conditions: 

a y u = 0{c- 4 ), A y u a = 0(c- 2 ), (28) 

where O y = r ) mn d m d n = d 2 + A y and A y = j eX d e d\ are the flat-space d'Alembert and Laplace operators in local 
{y m } coordinates. We shall use the term "harmonic metric tensor" to describe the metric tensor ([22"j) -([24" |) . expressed 
in terms of the harmonic potentials u and u a , given by Eqs. (1251) and (|26|) . which satisfy Eqs. (|2T[) and (|28|>. 



B. The form of the functions of the harmonic coordinate transformations 

As discussed in Sec. IIII A[ Eqs. (fT9]) and (|2"Tj) provide valuable constraints on the form of the metric tensor in a 
moving frame. As a matter of fact, these equations provide two additional restrictions on K, and Q a . It follows from 



4 An interesting consequence of the solution for the potentials u and u a is the fact that this choice ensures that the induced inertial 
space-time is flat, and thus, the Ricci tensor of this space-time vanishes, R m n{v) = 0- This can be verified by a direct calculation of the 
Ricci tensor R m n with the metric T] m n given by Eqs. 122B 1241 and also by applying Eqs. 11281 : 

R 00 = n y {c- 2 r,$ + C - 4 (r ? l 4 1 - i^ 1 ) 2 ) +0(c~ 6 )} = ~Dy{u + 0( C " 4 )} = 0{c~% 

R 0a = A s {c" 3 ^ +0( C " 5 )} = -7c,*^A s {« A + 0( C - 2 )} = 0( C - 5 ), 

R a p = A y {c- 2 V W +0(c" 4 )} = 7q/3 ^A !/ {« + 0( C - 2 )} =0(c- 4 ). 

This observation regarding the potentials u and u a will be important in the case of gravitational dynamics of iV-body systems, which 
will be studied in a subsequent paper. 
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Eqs. (|T9)) . (J2TJ) and the form of the metric tensor r) mn given by Eqs. (fTTj) (|T3)) that these two functions must also 
satisfy two first order partial differential equations: 



^ + «ba = 0( C - 4 ), (29) 



c^c^ +7qA ^ +7/3A 9^ +27q H^ + ^oJ = 0(c" )• (30) 

We now explore the alternative form of the harmonic gauge given by Eq. (|15|) . Substituting the coordinate transfor- 
mations (H])-© into Eq. (JTSJ), we can see that the harmonic gauge conditions restrict the coordinate transformation 
functions K, C and Q a such that they must satisfy the following set of second order partial differential equations: 

c 2lA +7 eA^i_ = 0(c - 2) (32) 

The general solution to these elliptic-type equations for the functions JC,C and Q a in Eqs. (j3"Tj) - (|3"3")l consist of two 
parts: a fundamental solution of the homogeneous Laplace equation and a particular solution of the inhomogeneous 
Poisson equation (except for Eq. pip , which is homogeneous). These solutions can be written in the form of a Taylor 
series expansion in terms of irreducible Cartesian tensors, which are symmetric and trace-free (STF) fl8j . 

The two sets of partial differential equations for JC, C and Q a given by Eqs. (|29[) -(|30" 1) and (|31|) (|33|) can be used to 
determine the general structure of these functions. 



1. Determining the structure of K, 



The general solution to Eq. (j3"Tj) with regular behavior on the world-line (i.e., omitting terms divergent when |y| — > 
or solutions not differentiable at |y| = 0) can be given in the following form: 

K{y) = K Q + K Gu y» + 8K + G{c- i ), where 5k = ^ -U Mfc (y V 1 '"^ + 0{c~% (34) 

fc>2 

with Kui...u k (y°) being STF tensors [l8j |, which depend only on the timelike coordinate y°. Substituting this form of 
the function K, into equation (|2U1) . we the find solutions for kq^ and K^...^: 

k „ = -cv „ + 0(c~ 4 ), = G(c- 4 ), k > 2. (35) 

As a result, the function /C that satisfies the harmonic gauge conditions is determined to be 

JC(y) = K -c(v 0u y' i ) + O(c- 4 ). (36) 

This expression completely fixes the spatial dependence of the function /C, but still has an unknown dependence on 
the timelike coordinate via the function Ko(y )• 



2. Determining the structure of Q a 



The general solution for the function Q a that satisfies Eq. (|33|) may be presented as a sum of a solution of the 
inhomogeneous Poisson equation and a solution of the homogeneous Laplace equation. Furthermore, the part of that 
solution with regular behavior in the vicinity of the world-line may be given in the following form: 

Q a (y) = qS + + + SC + o{ c -% (37) 

where can be determined directly from Eq. p3p and the function 5t; a satisfies the Laplace equation 



A y SC = 0(c~ 2 ). 



(38) 
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We can see that Eq. ((33]) can be used to determine (/(w, but would leave the other terms in the equation unspecified. 
To determine these terms, we use Eq. (f3"U)) together with Eq. (j2"9"|) . and get: 

dQ x dQ x /dlC \ 

v OaVop+1*\^+7e\^ + 2l a p[g-o + 2V 0e v e ) = 0(c- 2 ). (39) 



dy a ILip \dy° 

Using the intermediate solution (|36| for the function JC in Eq. (f39|) . we obtain the following equation for Q a : 



dQ x dQ x fdnQ \ 

«Oa"o/3 +7aA^j +7/3A^ + 2 7q/3 (^ + i WOe «o -a 0e y £ ) = 0{cr 2 ). (40) 

A trial solution to Eq. ([4"0"| may be given in the following general from: 

Q a = <Z Q + c^v^y* + c 2 v y y a + c 3 a°y e y< + c 4 a 0e y<y a + c 5 (|^ + ±v ^ )y a + y^ + S^(y), (41) 

where Qq and the antisymmetric matrix Wq 6 = — loq 01 are functions of the timelike coordinate y , c\, C5 are constants; 
and 8^(y), given by Eq. (|38|) . is at least of third order in the spatial coordinates y^, namely 5^(y) oc C(|y M | 3 ). Direct 
substitution of Eq. ([4"Tj) into Eq. (|4U|) results in the following unique solution for these coefficients: 

11 

ci = -~, c 2 = 0, c 3 = --, c 4 = 1, c 5 = -l. (42) 

As a result, the function Q a has the following structure 

Q a (y) = g Q -(^oX + < e + 7 Qe (^ + i^o A ))^ + % £ (yV-h Q W)+<5r(2/), (43) 

where q$ and Wq e are yet to be determined. 

By substituting (|43|) into (|40)) , we see that the function S£ a (y) in Eq. (|43p must satisfy the equation: 

+ d' 3 8C = 0{c~ 2 ). (44) 

We keep in mind that the function 8£ a (y) must also satisfy Eq. (f3"5)) . The solution to the partial differential equation 
with regular behavior on the world-line (i.e., when |y| — > 0) can be given in powers of y 11 as 

^ Q (y) = E l 5 ^, 1 ..., k (y°)y tll -y tlk + o(\ y »\ K ) + o( c - 2 ), (45) 



fc>3 



where S^q^ ^{y ) being STF tensors that depend only on timelike coordinate. Using the solution (|45j) in Eq. f|44[) . 
we can see that S^q^ „ fc is also antisymmetric with respect to the index a and any of the spatial indices 
Combination of these two conditions suggests that <5£oui Mfc = for all k > 3, thus 

<5r(y) - 0. (46) 
Therefore, application of the harmonic gauge conditions leads to the following form of the function Q a : 

Q a (y) = go Q -(H Q «o + ^o e + 7 Q lf^ + iW^ (47) 

where q$ , Wq £ and kq are yet to be determined. 



3. Determining the structure of L 

We now turn our attention to the second gauge condition on the temporal coordinate transformation, Eq. (|32p . 
Using the intermediate solution (|36[) for the function /C, we obtain the following equation for C: 

- - 2 fp + ^ = ^oA + - ^ (gj + + 0(c-'). (48) 
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The general solution of Eq. (|48|) for C may be presented as a sum of a solution SC for the inhomogeneous Poisson 
equation and a solution SCq of the homogeneous Laplace equation. A trial solution of the inhomogeneous equation 
to this equation, SC, is sought in the following form: 

SC = ck^aDiy^) + ck 2 {a Qe y%y v y v ) - fc 3C 2 ^(|^ + \v Qe vl)(y v y v ) + 0( C - 2 ), (49) 
where k\, k 2 , k 3 are some constants. Direct substitution of (|4^|) into (J3HJ) yields the following values for the coefficients: 

fa = ^, k 2 = ^, k 3 = i. (50) 
As a result, the solution for SC that satisfies the harmonic gauge conditions has the following form: 

SC(y) = I c ( Uoe ^)(^) + ^ C (W)^^ (51) 

The solution for the homogeneous equation (|48j) with regular behavior on the world-line (i.e., when |y| — > 0) may 
be presented as follows: 

My) = to(y°) + My ) y x + ! W2/ ) 2/ V + «(i0. (52) 

where £oa^ is an STF tensor of second rank and S£ is a function formed from similar STF tensors of higher order: 

S£ (V) = E ^W^O/V 1 "'^ + 0(\y»\ K ). (53) 

fc>3 

Finally, the general solution of Eq. (|48p may be presented as a sum of the special solution SC of the inhomogeneous 
equation and the solution Co of the homogeneous equation A y C = 0. Therefore, the general solution for the gauge 
equations for the function C(y) = Cq + SC has the following form: 

C(y) = 4+4a2/ A + ^oa m J/V + «(j/) + 

+ H v oA)(y,y v ) + ^cM)^y U ) ~ + \v 0e vi){y v y") + Q{c- 2 ). (54) 

We have determined the structure of the transformation functions /C, C, and Q a , which is imposed by the harmonic 
gauge. Specifically, the harmonic structure for K, is given by Eq. (|36)) . the function Q a was determined to have the 
structure given by Eq. (|48p . and the structure for C is given by Eq. f|54|) . Note that the harmonic gauge conditions 
allow us to reconstruct the structure of the functions only with respect to spatial coordinates The time-dependent 
functions Ko, (j^ , Wo % ^o,^ox^o\^, an d SIq^...^,. cannot be determined from the gauge conditions alone. We need 
to apply another set of conditions that would dynamically define the proper reference frame of a arbitrarily moving 
observer, thereby fixing these time-dependent functions. This procedure will be discussed in the following section. 



IV. DYNAMICAL CONDITIONS FOR A PROPER REFERENCE FRAME 



An accelerating observer that remains at rest with respect to an accelerating frame does so because of the balance 
between an external (physical) force that causes the observer to accelerate and the fictitious frame-reaction force that 
exists due to the choice of accelerating coordinates. The effects of both forces can be modeled in the form of an 
appropriately chosen effective metric rf^ n in the accelerating frame. In the case of complete balance, there will be no 
net force acting on the observer with respect to this effective metric, allowing it to be at rest in what we shall call its 
proper reference frame. The motion of the observer in this frame, then, will resemble a free fall that follows a geodesic 
with respect to the metric rf^ n . Thus, the observer's ordinary relativistic linear three- momentum, calculated in the 
accelerating frame and with respect to rf^ n , should be conserved. We can explore these conditions by writing down 
the Lagrangian of a test particle that represents the observer. By imposing further gauge conditions on the metric 
Vmni we nn( i that it is possible to eliminate all the remaining unknown components of /C, C and Q a . 
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A. Representing an external force using a fictitious harmonic metric 

We imagine an observer that remains at rest in the accelerating frame {y m }- We assume that this observer's 
acceleration b a is due to an external universal scalar potential that, in the vicinity of the observer's world-line, admits 
the following representation: 

<p(y) = <A) - e y e ), (55) 

where ipo — <fo(y ) is the external background potential on the observer's world-line and b a (y°) = —d a ip is the 
external acceleration acting on the observer. 

We incorporate the external potential into an effective metric, expressed in terms of y m . This allows us to use a 
single metric to capture all physical processes occurring near the world-line: the external potential, <p(y), that exerts 
the force accelerating the observer, and the frame potentials u and u a that generate a fictitious frame-reaction force 
that balances the force due to the external potential, keeping the observer at rest with respect to the accelerating 
frame. We shall denote this effective metric rf^ n . At the Newtonian order (or up to 0(c~ 4 ) in the t]qq component of 
the metric), we can write r)^ n by modifying the metric (|22 p -(|24 p as follows: 

ds 2 = (l - A(«(y) + p(y)) + 0( c - 4 )) (dy°) 2 + O(c- 3 )dy dy x + ( 7eA + 0( C - 2 ))dy*dy\ (56) 

In this metric, the potential ip is responsible for the external force and u is the inertial frame-reaction potential 
characterizing the accelerating reference frame, as introduced in the previous section and yet unknown. If the forces 
produced by the two potentials are equal to each other then a particle subject to the potential ip should not be 
accelerating with respect to its proper accelerating reference frame. This is the basic idea behind the method of 
formulating the dynamical conditions for a proper reference frame of an accelerated observer that we present below. 

We can see that the line element Eq. ([55)) is valid only at the Newtonian level. To extend this metric to the 
post-Newtonian level, first we write the acceleration as a sum of a Newtonian and post-Newtonian terms: 

b a (y°) = b^ a (y°) + c- 2 b [2]a (y°) + 0{c- 4 ). (57) 

Furthermore, we impose the isotropic harmonic gauge condition on the resulting metric: 

2cd oV [ ^ S + d^ eS = 0(c- 2 ), (58) 

[2]eff . [2]cff _2\ ,r n \ 

V a p + 7ctf%) = 0{c ). (59) 

This metric represents the combined contributions of the external force (introduced via b a ) and the inertial or frame- 
reaction force (introduced via u and u a ) that affect the motion of our observer in the accelerating frame. 

With the help of these equations we can reconstruct the general structure of the metric that corresponds to the 
acceleration (|57|) . It follows from Eq. (15811 that the acceleration-induced contribution to the mixed- index components 



of the harmonic and non-rotating metric rj^^ may be represented as 

vff = -^ a x(u\y)~\voy X + 0{y 2 )) +0( C - 2 ), (60) 

where tpa = cdoifQ. As we will be interested in the values of these potentials and their first spatial derivatives on the 
observer's world-line, we do not need an explicit form of the terms denoted by 0(y 2 ). Furthermore, Eq. (|59[) suggests 

that the spatial component of the metric tensor rj^ S in harmonic coordinates has the form: 

V [ T = 2 la0 (u{y) + V {y)) +0( C - 2 ). (61) 

The conditions lf6"0|) - (f6"Tj) allow us to extend the metric ([56"]) beyond the Newtonian level. In the fictitious metric 
rf^ n that is introduced to represent the combined effects of an external force and an accelerating frame, we must 
combine the contributions of the external and the frame-reaction forces: 

% e oG/) = l-^(u(y) + V (y))+^((u( y ) + v >( y )) 2 + 0(y 2 ))+0(c-% (62) 
r>fM = -lax^(u x (v) - ^o?/ A + 0(y 2 )) + 0(c-% (63) 

nfM = + («(») + <p(y)) + o^). (64) 



This metric represents accurately, at the post-Galilean level, the effects of an externally induced acceleration b a on 
an observer in an accelerating frame characterized by u and u a , while also satisfying the harmonic gauge. 
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B. Equation of motion of an accelerating observer 



The metric tensor r/^ n given by Eqs. (|62]) -([64 f allows us to study the dynamics of an observer or test particle that 
moves in response to the presence of the external force. The test particle Lagrangian L e g- that corresponds to this 
system can be obtained directly from the metric rf^n an d written as (l9j : 



L eS = -mc z — = -m<f\r), 



dy° V ,mn dy° dy° 



-mc 2 



{l + c- 2 ^v e v e -uj+ c- 4 (\u 2 - Hv^f - 4v e u ( + ^v e v e u + 0(y 2 )^j + 0(c~ 6 )}, (65) 



where u — u(y) + ip(y) + 0(c~ A ) and u a = u a {y) — hfioy a + C( c_2 ) are the combined scalar and vector potentials 
that consist of the inertial potentials given by Eqs. (f25|) . (|26[) representing the inertial frame-reaction potential, and 
the potentials Eq. (|55p combined with Eq. ([57]). representing the external force. 
The Lagrangian (p5|) leads to the following equation of motion 5 : 



d 

'dy~° 



— {v a (l + c- 2 (3u - \v,v e ) + G(c- 4 )) } 



-d a u{l-c- 2 (%v e v e + u)} + ^cd u a + ^v e (d e u a - d a u e ^j +lo{y)+0{c~ 4 ). (66) 



The condition that the test particle is to remain at rest with respect to the accelerating frame, then, amounts to 
demanding that its ordinary linear momentum be conserved, i.e., that its total time derivative is to remain zero. In 
other words, we require that there will be no forces acting on the observer in its proper reference frame or the right- 
hand side of (|6"6"|) vanishes on its world-line. Wc assume that the observer is located at the spatial origin, y a = 0. This 
leads to the following equation, constructed from the right-hand side of (|66l) that is valid on the observer's world-line: 

{ - d a u{l - c- 2 (lv t v € + uj } + ^cd Q u a + ^v e (>f t Q - d a u^ + ^0{y) + ©(c" 4 )} | _^ = 0. (67) 

We choose the coordinate system {y m } such that along the observer's world-line the potentials u and u a and their first 
spatial derivatives vanish, so that the metric rf^ n reduces to the Minkowski metric along the world-line. Therefore, 
we require that the following relations involving the frame-reaction potentials u and u a hold along the world-line: 

lim u(y) = -<p + G(c- 4 ), lim d p u(y) =b p + ©(c" 4 ), (68) 

lim u a (y) = 0(c~ 2 ), lim d p u a {y) = \5$ip Q + 0{c~ 2 ). (69) 

|y|->o |y|->o 

As we shall see in the next section, these conditions yield the equations needed to fix the time-like coordinate on the 
observer's world-line and to determine the explicit form of the coordinate transformation functions /C, C and Q a . 



C. Application of the dynamical conditions 



Imposing the conditions (j68|) - (|69|) on the potentials u and u a , which are given by Eqs. (]25 | - (j26|) . results in the 
following set of partial differential equations set on the observer's world-line: 

9kq -i , 1 f dC dQ e , /dn \ 2 /dn a , f \ 2 ~| _ 4 . 

<=o + <A) = - - - 2 vo,v - - j_ +CTOe _ + - 2 { w ) - [g-o + 2 v ev ) j = 0(c ), (70) 



5 It can be shown that Eq. H66H is equivalent to the geodesic equation written with respect to the metric rf^ n for the combined system of 
external and frame-reaction forces. However, Eq. I|66| l has the advantage that it allows to separate relativistic quantities and to study 
the motion of the system in a more straightforward way. Introducing this equation through the Lagrangian H65H offers us the opportunity 
to identify unambiguously dynamical quantities, most notably the canonical and ordinary momenta. Indeed we note that the canonical 
momentum is given by p" an = 8L e ff/dv a = v a (l + c~ 2 (3u — -^v e v eS j + 0(c~ 4 )) — 4c~ 2 fi a . In direct analogy with electromagnetism 
we note that the last term in this equation is an inerto-magnetic term; in contrast, the first group of terms corresponds to the test 
particle's ordinary (mechanical) momentum: p a = v a (l + c -2 (3ii — |« e i> e ) + 0(c~ 4 )) . 
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d u\ y=o -bp = a o0 -b f3 -^{b p + dy p dy0 + cv 0e Qy p dyQ + aop + v 0e v ) j = 0(c ), (71) 

1; iac 0G» 9Q A a 9« , n ,_ 2 , 



ltf=0 — 4 



1 d 2 C d 2 Q x d 2 Q x 

dpujy^-^ipo = -ho^o - ki- QycQyf) +claX dyVdyP +V0X dy a dyP + v ° aa °^ = ( 73 ) 



Eqs. (|70 p - (|73p may be used to determine uniquely the form of the coordinate transformation functions /C, C and 
Q a . Indeed, from the first two equations above, (|70|) and (|7T|) . we immediately have: 

§«^§ = 0{c-% (74) 
a° = b^ a +0(c- 4 ). (75) 
Using Eq. (|74p in Eq. ([55)1 , we can determine uniquely the function /C: 

Substituting this expression into Eq. (|47p we can determine the function Q a : 

Q a (y) = ^-(5<«o + ^o e + 7 Q Vo)ye + ao e (yV-57 ae yAy A )+0( c - 2 ). (77) 

Finally, the general solution for the function C given by Eq. ([54]) now takes the following form: 

C{y) = to + lox V X + |W 2/ V + gc((«oe<»o) - V^o) (&,!/") + 1 7 5 c(a 0e 2; e )(y,2/ ,y ) + + 0(c- 2 ). (78) 

The next task is to find the remaining undetermined timc-dcpcndcnt functions present in Q a and C, as given by 
Eqs. (|77jl and (|78|) . To do this, we rewrite the remaining parts of Eqs. ([70 |) ~(f73 |l as a system of partial differential 
equations with respect to C, again set on the observer's world-line: 

d£ 9Q e x fdKo\ 2 fdn i e \ 2 _ 2 , 

w +CVOe w + ^w) -\w +M ) ( } ' (7) 
w +CTo w +aw (^ + ^) = - & " 1+c(c ~ 2) < (80) 



^ 9Q A | dQ a va dKp 
c dy e 7 ° A dy € dy° dy° 



r e ~irz+r e v x^- + c^--vZ-Z = 0{c~% (81) 



ld^ + Vox ^^ + claX M? + VOaaop = -h^m + o{c- 2 ). (82) 

The equations above are may be used to determine the remaining unknown time-dependent functions £q, £o\, (-oxp,, and 
also q$ and Wg e still present in the coordinate transformation functions. Thus, substituting the previously obtained 
solutions for K, and Q a , given by Eqs. (|76p and (|77p. in Eq. ([75]) leads to the following solution for £o : 

iio = -v eq e o- l^v^f + \{v 0e iQip + \ V l + O{c- 2 ). (83) 

Next, Eq. (|80|) results in the equation for ^ : 

|iS" = -b l2]a + H Ka ) " + ft - «o + 0(c~ 2 ). (84) 

From Eq. ||SIJ) we can determine £g : 

= -g tt - t ; 0eWo Qe + 2 t ; V +0( C - 2 ). (85) 
Eq. (|82p leads to the following solution for 1^: 

\lf = -ftfag - i^a" + | 7 ^(a 0e ^)+d; ^ + O(c- 2 ). (86) 



12 



The quantity £q is an STF tensor. The expression on the right-hand side must, therefore, be also symmetric. This 
can be achieved by choosing the anti-symmetric tensor Wq ^ appropriately. This can be done uniquely, resulting in 

^ = K^-^)+0(^), (87) 
which is the relativistic Thomas precession (2TJ]. Using this expression yields the following solution for 

\lf = -ttfog - «£o? + §7 Q V«o + 0{c~ 2 ). (88) 
Furthermore, we can use Eq. ([57)) in Eq. ([8"4")) . leading to the solution for £q in the form: 

\t% = -bW° + ±a%(v 0e v< ) + vgtpo - ogVo + 0(c- 2 ). (89) 
Eqs. ([85]) and ([89)) allow us to determine q$ . Specifically, rewriting Eq. ([89]) as 

ki° = - b m° + i(<(« 0e ^))' - « Q (« X) + (« >o)' - 2a^ + G(c- 2 ) (90) 
and formally integrating it with respect to y°, we obtain another expression for £q : 

W - l«."(».,«5) + ** - / + »."(».,«5) + iV + o (c->). (91) 

Eqs. (|85p and ([9T|) can now be solved with respect to <?q : 

<Zo = ~ l 2 VoM) ' v e"o £ + <</>o + / + «?(« 0e 6 [0]6 ) + 2&[0l >o) + 0(c- 2 ), (92) 

where we used Eq. (j75[) for Oq. The first term in this expression is the Lorentzian factor, while the other terms 
explicitly depend on the observer's acceleration and the value of the external potential tp = (fio(y°) on its world-line. 

Finally, differentiating Eq. (|85[) with respect to time and subtracting the result from Eq. (|89[) (or just simply 
differentiating Eq. (|92[) with respect to time), we obtain the following equation on : 

qg = b^ a - + co^aoe + Mg<Po + + 0( C - 2 ). (93) 

Therefore, with the knowledge of the external acceleration b a and background potential tpo, we can use Eqs. (|92p and 
([93)) to completely determine the function qfi . 

The true position of a test particle includes terms to all orders, not just the first-order (Galilean) term. This led 
us to introduce the vector, Xq(u°), defined by Eq. (j?)). Combining this definition with Eqs. ([75)) and ([53")h we can 
now write the magnitude of the frame-reaction force (acting on the unit mass) written in the local coordinates of the 
accelerated observer as it relates to the measured acceleration b a : 

*?(V ) - «o + c~ 2 9o + 0(c- 4 ) = b® a + ±{b^ a - (fttfvg + u^b® + 3afa„ + ttfft,} + 0(c" 4 ) 

= &e{7 Qe -^(^o^+< e )} + ^{3a^o + «oVo}+0(c- 4 ). (94) 

The equation of motion ([M)) establishes the correspondence between & Q , the externally- induced acceleration of the 
observer, and the fictitious frame-reaction acceleration Xq (y°) that is needed to keep the observer at rest in its proper 
reference frame. This frame-reaction force balances the effect of the external inertial force acting on the observer. 

D. Summary of results for the direct transformations 

We sought general post-Galilean transformations between the Minkowski frame {x k } and the dynamically non- 
rotating coordinates {y k } of a proper reference frame of an accelerated observer. We did that representing such a 
coordinate transformation in the most general form: 



x° = 2 / + C - 2 /C(y , 2 / e )+ C - 4 £(y , 2 / e ) + O(c- 6 ), 



(95) 
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x a = y a + z«(y ) + c- 2 Q a (y°,y<) + O(c- 4 ). (96) 

To determine the unknown functions /C, C, and Q a , we used the following approach: i) we imposed the harmonic 
gauge conditions on the accelerated Minkowski metric, the metric tensor in the local accelerating frame; ii) we ensured 
that the accelerating frame is non-rotating and the chosen coordinates are spatially isotropic; iii) we introduced an 
accelerating observer at rest with respect to its proper accelerating frame; iv) we cast the combination of all forces 
acting on the accelerating observer in the form of a fictitious metric that includes the external and frame-reaction 
potentials; v) we imposed the harmonic gauge on the fictitious metric; vi) we required that a co-moving test particle's 
ordinary three-dimensional linear momentum be conserved on the world-line occupied by the accelerating frame. 

Together, these conditions were sufficient to determine K., C and Q a unambiguously: 

= / (<Pa- h v oA)dy'° - c(v ^) + o( c - 4 ), (97) 

£(y) = -/ (l'0e<?0 + |( U Oe«o) 2 ~ K^Oe^o)^ ~ 3Vo)<V° ~ c(qOe + «0 w 0eA ~ 2v 0e ^ )y e - 
Jyo 

-%c(a 0e v ox + a ox v 0e - "fexa 0tl vg + ±7^0) 2/ V + ^c(a 0e y e )(y v y u ) + 5£{y) + 0{cT 2 ), (98) 

Q a (y) = <Zo -(^o«o + < £ + 7 Q Vo)^ + «o £ (y Q y £ -h Q W)+0(c- 2 ), (99) 

with the anti-symmetric relativistic precession matrix ujq a given by Eq. (j87[) . and the post- Newtonian component of 
the spatial coordinate in the local frame, q$ , given by Eq. (|93l) . 

Substituting these solutions for the functions /C, C, and Q a into the expressions for the inertial frame-reaction 
potentials u and u a given by Eqs. (|2"5j) -([2l) |) . we find the following form for these potentials: 

u(y) = {a^)- Vo + ^[\(iy'y x -Y X y^)a 0t a ox -^{a^){y (100) 

u "(y) = -^{Sy a y £ -i a V)a 0e + ^ y a ~d a i- c se + o(c- 2 ), (101) 

where a a denotes the frame-reaction acceleration which is equal to the measured external acceleration b a given by 
Eq. (|57j) or a a (y°) = b a {y°) = 6[°) Q + c~ 2 b^ a +0(c" 4 ). Substituting these expressions for the inertial potentials into 
Eqs. ([2"2" j) - (j2"4")) leads to the following form of the accelerated Minkowski metric of the arbitrarily moving observer: 



+ ^(aocV'KVuV*) ~ kMVrV*) + + <?( c ~ 6 )- (102) 

Vo a (y) = 7aA^{^(3yV~7 A %y' i )ao e -^oy A } + ^a Q ^ + 0( c - 5 ), (103) 

ri a p{y) = 7a/3+7c t /3^{(a 0e 2/ e )-^o}+C , (c" 4 ). (104) 

All terms in this metric are determined except for the function S£, which remains unknown. We note that the potentials 
Eqs. (|100[) - (|101[) depend on the partial derivatives of 6£. The same partial derivatives appear in the temporal and 
mixed components of the metric (|102[) - (|103[) . The presence of these terms in the metric amounts to adding a full time 
derivative to the Lagrangian that describes the system of the moving observer. Indeed, separating in the Lagrangian 
constructed from Eqs. (|102|) - (|104j) the terms that depend on 5£, we have: 

As a result, the remainder of the gauge transformation depending on S£ will not change the dynamics in the system 
and, thus, it can be omitted. After some re-arranging, the frame-reaction potentials u and u a take the form: 

u(y) = {a^)-Vo + ^{\(^a Qt a QX ~ ltX a Q ^ + \ ltX ^y e y X ~ ^ (106) 
u a (y) = -^(3y a y*- 7 a V)a 0e + &oy a + C(c- 2 ). (107) 
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We can now present the metric of the moving the observer in the following final form: 

Voo(y) = 1 - ^{(aey e ) - <A)} + Jf{(( a 0e2/) - Vvf - 

- h( 3a Oe a o\ ~ 7 e Aa OM a o + gTeAPoJj/V + To (a 0e y £ ) (V^) } + C( c ~ 6 )> 

Vaf)(y) = la/3 +l a f}^[(a 0e y e ) - ^0} + C(c~ 4 ). 
The coordinate transformations that place the observer into this reference frame are given below: 

x° = y° + c- 2 {[ ("^ -I(« 0£ + c -% e )(^ + c - 2 ^) + c - 2 (I^ + I( UOei; e )¥3o _ ^ VOeV ^\ dy >o + 
-c(vl + c- 2 g ) ( 7 ,a - c- 2 (uj QcX + 2 7£A ^ ))y A } - 

~ C ~ 4 {l C ( a 0e V 0\ + a 0\ v 0e ~ 7eA<iO/X + §7eA¥>o)yV ~ ToC^y'Xy^)} + C(c~ 6 ), (HI) 

x a = y a + z« + c- 2 { q %- (K v o + < e + 7 Q Vo)y e + «oe(y a 2/ e -h Q£ yAy A )} + 0(c- 4 ). (112) 

The presence of ipo m the metric tensor (|108p - (|110[) and coordinate transformations i|112p is quite interesting. 

It shows that in the case of a time- varying background potential <po(y°) an d no external acceleration, the metric tensor 
of the corresponding space-time differs from the Minkowski metric. Furthermore, the new space-like coordinates are 
scaled by tpo and the time- like ones are stretched by both ifo and (f>Q. The result may be intuitive, but was not 
available previously. One can verify that in the case of uniform constant velocity motion (tig = 0) and in the absence 
of the external background potential (ipo = 0), the metric given by Eqs. (|108p - (|110[) reduces to the Minkowski metric, 
Vmn = Imn- Also, according to Eq. ([52"j) and setting ip = 0, the function becomes <jg = — ^vfi (v 0e v^ + 0(c~ 2 ), 
and the transformations Eqs. (|lll|) - ()112j) above reduce to the Lorentz transformations. 

The expressions (|108|) - (|110p represent the harmonic metric tensor in the local coordinates of the accelerating ref- 
erence frame. This metric and the transformations (jllip (|112l) are new and extend previous formulations obtained 
with different methods. We were able to derive for the first time an explicit form of the metric tensor corresponding 
to the space-time of an accelerated observer under harmonic gauge conditions and corresponding coordinate transfor- 
mations. These results may be verified in laboratory conditions, for instance, those involving high-energy accelerators 
or precision physical measurements. The formulation can be used to develop models for high precision experiments 
(for example, those discussed in (ll I2I [2lT[23| ) where one would needed to relate various observable quantities that arc 
critical for experimental success. However, for a complete description of these experiments we would need to establish 
the inverse coordinate transformations - the task that will be performed in the next section. 

V. INVERSE TRANSFORMATIONS 

In the preceding sections, we constructed an explicit form of the direct transformation between inertial and acceler- 
ating reference frames by applying the harmonic gauge and dynamical conditions on the metric. Given the Jacobian 
matrix it was most convenient to work with the covariant form of the metric tensor, which could be expressed 

in terms of the accelerating coordinates by trivial application of the tensor transformation rules (|10l) . 

The same logic suggests that if we were to work on the inverse transformation: that is, when it is the inverse 
Jacobian matrix dy m /dx n that is given in explicit form, it is more convenient to work with the contravariant form of 
the metric tensor, to which this Jacobian can be applied readily. This simple observation leads us to the idea that we 
can get the inverse transformations — i.e., from the accelerated to the inertial frame — by simply repeating the previous 
calculations, but with the contravariant form of the metric tensor instead of the covariant form. 

In this section, we show that this is indeed feasible, and accomplish something not usually found in the litera- 
ture: construction of a method that can be applied for both direct and inverse transformations between inertial and 
accelerating reference frames at the same time, in a self-consistent manner. 



(108) 
(109) 
(110) 
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A. General form of the post-Galilean coordinate transformations 

We write the inverse of the general post-Galilean transformations between the dynamically non- rotating 

coordinates of accelerating {y m } and those of inertial {x m } reference frames in the following form: 

? y° = x° + c- 2 /C(x ,x e ) + c~ 4 C(x lX e ) + O(c- 6 ), (113) 

y a = x a -z«{x ) + c- 2 Q a (x°,x e ) + O{c~ 4 ), (114) 

where Zq(x°) is the Galilean vector connecting the spatial origins of the frames, expressed as a function of global time, 
x . Our objective is to determine the functions fC, £ and Q a in explicit form. 

We can verify that, in order for Eqs. (fTT3| -(fTT4 |) to be inverse to Eqs. ([S])-©, the "hatted" functions {JC,C, Q a ) 
must relate to the original set of (IC,C, Q a ) via the following expressions: 

t{x) = -K,(x ,r e ) + O(c- 4 ), (115) 

t{x) = d K.(x ,r e )-K.(x ,r e )+d x K.(x ,r e )-Q x (x Q ,r e )-£( X °,r e )+O(c- 2 ), (116) 

Q a (x) = {v^/c)lC{x y)~Q a {x n y)+O{c- 2 ), (117) 

with r e = x e — Xq, where Xq(x°) is the post-Galilean vector between the origins of two non-rotating frames expressed 
as a function of the global time-like coordinate x° (as opposed to Eq. which is given in local time y°) defined as 

x% = 4 - c- 2 Q%t°,0) + 0(c- 4 ), (118) 

and also, Vq = ip , do = d/dx° and d\ = d/dx x , and x = (x°,x e ). 

The inverse of the Jacobian matrix ([T])-©, dy n /dx m , can be obtained directly from (|113p - (|114[) : 

8x° + dx" + dx 0+ { h dx» dxf* + dx^ + 1 h 1 ' 

W v& a d& s dy» 2 d& t 

W = -~ +C 9^ + ° (c } ' dx~^ d » + C ^ + ° {C ] - (120) 

We note that this Jacobian matrix is composed of the quantities dy m /dx k that are clearly functions of {x k }. 
Using the Jacobian matrix and standard tensor transformation rules, we can express the relationship between the 
contravariant Minkowski tensor and the contravariant metric of the accelerating frame in the form 

*™w = -hw^ x(3, » or ^ mn{y ( x)) = i^^ kt{x) - (121) 

We denote rj mn (x) = r] mn (y(x)) and, taking the somewhat unusual step of using the contravariant tensor transfor- 
mation rule (|121[) together with Eqs. (|119[) - (|120[) . we obtain explicit expressions for the contravariant components of 
the accelerated Minkowski metric ff nn {x), expressed as functions of the global coordinates {x k }: 

= ^ + ?{^o+7"^+7^fr} + 0(c- 4 ). (124) 

As in the case of the direct transformation, we impose the harmonic gauge condition on the metric, to help us 
establish explicit forms of the transformation functions t, C, and Q a . 



B. Imposing the harmonic gauge condition 

Similarly to the case of direct coordinate transformations, we will use the harmonic gauge conditions. In analogy 
with the derivation of Eqs. (|22|) -(f24 | . we derive the contravariant metric f) mn : 

77 00 (x) = l + \u{x) + ^u 2 {x)+0{c~ 6 ), (125) 



16 



fj 0a (x) = ^u a (x) + 0(c- 5 ), (126) 
ff^( x ) = 7 ^- 7 ^4u(x)+0(c- 4 ), (127) 
where the inertial scalar u(x) and vector u a (x) potentials are now expressed via the "hatted" transformation functions: 

«M = S5 + i-»-"S + Aw + fa? + Ha?) - (a? + '"»'" s ) } + ° (c >■ (128) 

One can verify that these two potentials satisfy the following continuity equation in global coordinates {x m }: 

(cd + v e Q d e )u + d e u e = 0(c- 2 ). (130) 

In addition, by a direct calculation with the help of Eqs. (|135[> — (|137|) and (|132p - (|133[) . one can verify that these 
potentials also satisfy the harmonic equations: 

U x ii = 0{c~% A x ii a = 0(c- 2 ), (131) 

where D x = j mn d m d n and = j eX d € d\ are the d'Alembertian and Laplacian, correspondingly, with respect to 
{x m }. As before, we refer to the metric tensor expressed in terms of u and u a as the harmonic metric tensor. 

Although the expressions for the scalar inertial potentials have different functional dependence on the transformation 
functions (i.e., (1C,£, Q a ) vs. (tt,£, Q a )), it is clear that the expressions for u(y) and u(x), given by Eqs. (|25|) and 
(|128|) . represent the same quantity that is being expressed in terms different coordinates: local {y m } and global {x m }, 
so that u(x) = u(y(x)). The same is true for the inertial vector potentials u a (y) and u a (x), given by Eqs. (|26[) and 
(|129p . for which the following relation holds u a (x) = u a (y(x)). 



C. The form of the functions of the harmonic coordinate transformations 

In addition to the constraints provided by Eqs. (|19|) and (|21[) on the form of the metric tensor in a moving frame, 
we can once again derive two additional equations on K, and Q a . From Eqs. (|19p and (f2Tj) and Eqs. ()122[) - (|124|) . we 
find that these two functions must also satisfy two first order partial differential equations: 



1 dK 

c dx £ 

^+7-^ + /^ + 2 7Q ,(— + I 7 ----^) = 0(c-»). (133) 

In analogy with the derivation of Eq. (| 15[) given in Sec. lIIIBl we now use the harmonic gauge given by the equation: 

U lV m = 0, (134) 

where D 7 denotes the covariant d'Alembertian with respect to the metric of the inertial frame 'y mn ( x )j acting on 
y m {x). Substituting the coordinate transformations (|113[) - (|114[) into this equation, we can see that the harmonic 
gauge conditions restrict the coordinate transformation functions IC, C and Q a only to those that satisfy the following 
set of second order partial differential equations: 

^ A aS? " ^>- < 135 > 
-< + ^' X ££x = o [c r% ,137) 

The two sets of partial differential equations on IC, C and Q a given by Eqs. (|132p ~ (|133|) and (|135p (|137|) can be used 
to determine the general structure of these functions. In strict analogy with Eqs. (|36|) . (|47|) and (|54|) . we can derive 
the following form for the functions K,,C and Q a : 

£{x) = k + c(v ^) + O(c- 4 ), (138) 



7 Q£ -^7-< = 0(cr% (132) 
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C(x) = to + lox r A + ±£ 0V rV + \c(v Qe al){r v r u ) - i c (a 0£ r e )( V) - 

^ c2 ^ (S ~ ^) + 5i{x) + o(c_2) ' (i39) 

Q"(x) = -^-(K^ + C + Tl^-M)^%(^47V)+0( C - 2 ), (140) 
where r a is defined by 

r a = x a ~ x%, x% = z% + c~ 2 q« + 0(c- 4 ), (141) 

and Xq (x°) is the vector that connects the origins of the two reference frames, complete to 0(c~ 4 ) and expressed as a 
function of the global time-like coordinate x°, given by Eq. (|118[) . We can see that k ,qQ,ui^ e , ^Oj^oAj^oA/* an d S£(x) 
cannot be determined from the gauge conditions alone. Similarly to the case of direct transformations, we need to 
apply another set of conditions that would fix the reference frame of a moving observer, which we discuss next. 



D. Finding the form of the coordinate transformation functions 

In the case of the metric tensor given by the expressions Eqs. ()125[) - (|127j) . the approach discussed in Sec . II V Al and 
HVBI that yielded the conditions Eqs. (|68| -(|69 ]l . similarly leads to the following set of equations: 

lim = -^ + C(c" 4 ), lim dpu(x) = bp +C(c" 4 ), (142) 

|r|— >0 l r |->0 

lim ii a (x) = 0(c- 2 ), \im dfsu a (x) = \61<i> + O(c- 2 ), (143) 

|r|->0 |r|->0 o P' 

where, similarly to Eq. ([57} , ipo(x ) = (po(x° + c~ 2 ko + 0(c~ 4 )) is the background potential at the observer's world-line 
(where we used y°(x a ) = x° + c~ 2 ko(x ) + 0(c~ 4 ), valid on the world-line) and the measured acceleration b a is related 
to the coordinate acceleration Xq , given by Eq. (|94)l . as b (x°) = (l — 2tpo/c 2 + 0(c~ 4 ))xq (y°(x )) , so that b relates 
to b a as 

b a (x°) = 6^x + C - 2 ^ + 0( C - 4 )){ 7 ae -l(i<^+< £ -7 Q Vo)} + ^ Vo + 0(c- 4 ). (144) 

Imposing the conditions Eqs. (|142[) - (|143[) on the potentials u and u a given by Eqs. (|125|) - (|127[) results in a set of 
partial differential equations set on the world-line of the accelerated observer. These equations can used to determine 
the coordinate transformation functions entering Eqs. (|138|) - ()140|) which were found to be: 

ko = - f (<p - |u e«o - c~ 2 v 0e q^jdx m + 0(c~ 4 ), (145) 
< = b^ a + 0{c- 4 ). (146) 
In addition, the use of Eqs. (|142j) - (|143[) yields the following solutions for the functions £o, (q, fp, and £q a : 

(147) 
(148) 
(149) 
(150) 

Using the same argument concerning the symmetry properties of £q ^ that led to Eq. (|87jl , we find the following unique 

x a/3 . 

choice for the anti-symmetric matrix lu representing the Thomas precession: 

Qf = -u,? = -i{v^ -v^)+O(c- 2 ). (151) 

Finally, Eqs. (|150|) and ([59")) allow us to determine the equation for q^ . Indeed, differentiating Eq. (|150[) with respect 
to time and subtracting the result from Eq. (|149[) . we obtain: 

% = b^ a + + agvftvto + vSfo + 0(c~ 2 ). (152) 





= -\{v 0t v c ) 2 - \{v 0e vl)(p -1 


-I0l + O{c~ 2 


-M 


= & [2]Q + ~ <fa)a$- 


f 0(c- 2 ), 




= % - f l, o(%«o) -<0o + 


0(c- 2 ), 






+ 0(c~ 2 ). 
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By combining (|146|) and (|152[) . wc obtain the equations of motion of the accelerated observer with respect to {x m }: 
i Q (*°) = a% + c- 2 C + 0(c- 4 ) = b® a + ±{ti 2 ^ + (« + a%v' )v 0e + vg0 o } + 0(c~ 4 ) 

= ^{^ + ^K«o +7 a£ ^<)} + ^«o0o + G(c- 4 ). (153) 

The equation of motion (| 153[) establish the correspondence between the external force acting on the test particle, 
b a , and the resulting acceleration of the moving frame Xq (x°) as measured by the inertial observer. 



E. Summary of results for the inverse transformation 

The coordinate transformations between the global coordinates of the inertial frame and the local coordinates 
introduced in the proper reference frame of an accelerating observer have the following form 

y Q = x° + c- 2 IC(x ,x e ) + c- 4 C(x 1 x e ) + O(c- 6 ), (154) 

y a = x a - z«(x ) + c- 2 Q a (x a ,x e ) + O{c- 4 ), (155) 

with the transformation functions 1C,C and Q a given by 

t{x) = - f (ipo - " c~\Ao)dx'° + c(v Q y) + 0(c- 4 ), (156) 

= I {- U v Oe v o) 2 ~ ^( v OeK)^o + ^l)dx'° + c(q ox - ^v ox {v 0e v^) -v ox ^o)r x ^ 

+ H v ne a ox + vox^oe + §7eA0o)r e r A - ±c(a Qe r e )(r u r v ) + S£(x) + G{c~ 2 ), (157) 
Q a (x) = -C- (H aw o + ^o e -7 Q ^o)r- £ -ao e (r- Q r e -i 7 Qe r A r A ) +0(c- 2 ), (158) 

where r a = x a — Xq is given by Eq. (|14ip , the anti-symmetric relativistic precession matrix u>q a is given by Eq. ()151[) , 
and the post-Newtonian component of the spatial coordinate origin of the local frame, is given by Eq. (|152[) . 

Substituting these solutions for the functions /C, C and Q a into the expressions for the inertial potentials u and u a 
given by Eqs. (|128p ~ (|129p . we find the following form for these potentials: 

u{x) = (a e r £ ) - <p + ^{^{a Qe v 0X + v Qe a 0X + a 0e a 0A + ±7^0)^ - ^(aoe^XV')} + 

Mcdo + v^d^Si+Oic- 4 ), (159) 

= -±{3r a r e -rV}a 0e + ^ r a +d a j- c 5£ + O(c- 2 ) : (160) 

where a a denotes the frame-reaction acceleration, which is equal to the external acceleration b given by Eq. (|144[) 
and measured in the global frame; i.e., a a (x°) = b a (x°) = b^ a + c~ 2 U 2 ^ a + 0(c~ 4 ). 

The same argument that allowed us to eliminate 81 in Sec. IIVDI works here, allowing us to omit 81. Therefore, 
we can now present the metric of the moving observer expressed in the global coordinates {x"}. Substituting the 
expressions for the inertial potentials given by (|159|) and (|160|) into Eqs. (|125j) - (|127[) . and omitting S£, we have: 

= l + |{(a e r e )-^ } + |{(( a0e r e )-^) 2 + 

+ h( h neV \ + v 0t a ox + a ae a ox + §7^0) r e r* - ^(a j- £ )(r M r' 1 )} + C>(<r 6 ), (161) 



f, 0a (x) = --{i(3r<V _ 7 «V0do e - & r°} + 0{cr% (162) 
f(*{x) = 7^ - 7 ^l{(a 0e r e ) - ^0} + 0(c- 4 ). (163) 



2 

c 

The coordinate transformations that put the observer in this reference frame are given by 



y° = x° + c 2 { h{v Qe + c 2 i 0e )(v t + c 2 q e ) + c 2 {i(vS) 2 +lM)^-^o})^'° + 
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+ C (l - C- 2 (±( V £) + ft))) (« 0e + c" 2 ^e)?- £ } + 

+ c~ 4 {^o £ aoA + «0e«0A + H^o)r e r A - ^(do^^^r^)} + 0(c- 6 ), (164) 
y" = r tt - C - 2 {(i^+^- 7 a ^ )r £ +ao e (r Q r e -i 7 Q V A r A )}+0( C - 4 ), (165) 

where r a = x a -x%{t), as defined by Eq. (flUjl . Eq. (|TT7)) suggests that = cft-(vg/c) J((f - \v Qe vl)dy'° + C(c" 2 ), 
which, after accounting for Eq. (|92[) . in the case of constant velocity motion and the absence of the background external 
potential tpo = 0, results in q — 0(c~ 2 ). Consequently, in this case the metric given by Eqs. (|161|) - (|163|) reduces 
to the Minkowski metric, r) mn = <y m ™, and the transformations above reduce to the Lorentz transformations between 
the coordinates of the accelerated (y°,y a ) and inertial (x°,x a ) reference frames: 

V° = (l +c- 2 i( % ^) - c- 4 !^^) 2 )! + (l - c- 2 !^^)^ 1 ^'' 6 ) + 0(c' 6 ), (166) 
y a = (^-c->o% e )r £ + 0(c- 4 ). (167) 

If we now write x[)(t) = VqC x , wc obtain expressions identical to Eqs. (|5])-([6|), demonstrating correspondence of 
our results to those established previously using different techniques. 

This set of results concludes our derivation of the coordinate transformations from local coordinates of an accelerated 
reference frame to the global coordinates introduced in the inertial frame. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we introduced a new approach to construct coordinate transformations between inertial and accel- 
erating reference frames. Our objective was to establish properties of a frame associated with a given accelerating 
world-line and to find the metric tensor corresponding to this frame; wc achieved this goal by imposing a set of clearly 
defined coordinate and physical conditions. Specifically, we combined the use of the harmonic gauge and physical 
insight based on the dynamical properties of the proper reference frame of an accelerating observer. This way, we were 
able to find a unique non-rotating, accelerating coordinate system that corresponds to a given accelerating world-line. 
We derived an explicit form of the coordinate transformation between the inertial and accelerating reference frames 
to the first post-Galilean order. The transformations remain valid so long as v <C c. In addition to the direct trans- 
formation, which allows us to express the coordinates {a;" 1 } of the inertial reference frame in terms of the accelerating 
coordinates {y m }, we also developed the inverse coordinate transformations in explicit form. 

The metric associated with t he p roper reference frame of an accelerated observer has been studied by many re- 
searchers (see for instance, [iMa, |24| ) . It was found that when only the acceleration Oq of the reference frame is taken 
into account, the corresponding metric in local non-rotating coordinates {y' m } = {y'°,y' a ) takes the form pi [25|: 

% s o(y') = (i + vl(y') = o, v%{v') = 7«/»- (168) 

To establish the correspondence of the newly found metric r] mn (y) of an accelerated observer in harmonic coordinates 
of its proper reference frame Eqs. (|108[) - (|110p to the metric Vmniv') gi ven m the form of Eq. (|168|) we can use the 
same approach employing the functional /C£Q-parameterization. To do this, first of all, we note that in the absence 



of the external potential (ipo = 0) the metric Eqs. (|108|) (|110[) reduces to: 

Voo(y) = 1- jr(a £ y e ) + Jfl^AaoX-ao^oA^V + (169) 

Voa(y) = 7aA^{^(3yV-7 A W>o e }+0(c- 5 ), (170) 
2 

Vapiv) = 7a/3 +7c/3 — (a 0e y £ ) + C(c -4 ). (171) 



Then, in analogy to Eqs. ([I])-© we can introduce coordinate transformations y rn = f ,n (y' k ), which would depend 
on a set of to-be-determined functions of IC\C, and Q' a . Matching the two metric tensors via the usual tensor 
transformation rule rj mn (y) = (dy' k /dy m )(dy n / dy n )ri^ l (y') yields unique solutions for these functions. As a result, we 
can obtain the following coordinate transformation that transforms the well-known metric of an accelerated observer 
in coordinates (|168p to the new harmonic accelerated metric given by Eqs. (|169[) - p7ip : 

y'° = V° + c- 4 {c±(d 0t y*)(y x y>)}, y' a = y a + C - 2 {a 0e (y a y* - ± 7 Q W)}- (172) 
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Note that the inverse transformations are obtained simply by replacing (y'°,y' a ) with (y°, y a ) and the acceleration <Xq 
with — <2g . This establishes correspondence between our results reported in [3U5l.l24| that were previously derived under 
different gauge conditions and using different techniques for constricting the proper reference frame of an accelerated 
observer. The resulting metric (|169p - (|17ip can be used to study physical phenomena in this frame. 

The new results reported in this paper find a good correspondence with the results previously obtained by other 
researchers who also used the harmonic gauge conditions. Thus, the coordinate transformations that we derived in 
the context of a flat space-time find their exact correspondence with the results established in the presence of gravity 
for both direct and inverse coordinate transformations given in 043 and @, Ill"l - fl3| correspondingly. However, our 
approach allows one to consistently and within the same framework develop both direct and inverse transformations, 
corresponding equations of motion, and explicit forms of the metric tensors in the various reference frames involved. 
The difficulty of this task was mentioned in [To| when the post-Newtonian motion of a gravitational A^-body system 
was considered; our proposed formulation successfully resolves this important issue. As an added benefit, the new 
approach provides one with a good justification to eliminate the functions 5k, 51; and 51 yielding a complete form for 
the transformation functions K,, £ and Q a involved in the transformations (as well as their "hatted" counterparts). 

The significance of our result is that for the first time, a formalism for the coordinate transformation between inertial 
and accelerating reference frames is provided, presenting both the direct and inverse transformations in explicit form. 
By combining inverse and direct transformations, the transformation rules between arbitrary accelerating frames 
can be obtained. Furthermore, it is possible to combine direct (or inverse) transformations, and obtain another 
transformation that can be represented by our formalism, as shown explicitly in |17| . This leads to an approximate 
finite group structure that extends the Poincare group of transformations to accelerating reference frames. 

The results obtained in this paper provide a clear framework to describe observables, field transformations, and 
corresponding equations of motion that are needed to describe modern-day high-precision experiments. Specifically, 
the new approach is designed to facilitate the analysis of relativistic phenomena where effects of acceleration may be 
significant. We should note that the approach we presented can be used in an iterative manner: if greater accuracy 
is desired, the coordinate transformations (HJ)-© can be expanded to include higher-order terms. Furthermore, the 
same approach relying on the functional /C£Q-parameterization may be successfully applied to the case of describing 
the gravitational dynamics of an astronomical A^-body system and dynamically-rotating reference frames. This work 
has begun and the results, when available, will be reported elsewhere. 
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